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ABSTRACT 

Semi empirical methods have shown considerable promise in aiding in the calculation of many prop- 
erties of materials [1,2]. Materials used in engineering applications have defects that occur for 
various reasons including processing [3]. In this work we present the first application of the BFS 
(Bozzolo Ferrante and Smith) method for alloys [1] to describe some aspects of microstructure due 
to processing for the Cu-Au system (Cu-Au, CuAu 3 , and Cu 3 Au). We use finite temperature Monte 
Carlo calculations, in order to show the influence of ’heat treatment’ in the low-temperature phase 
of the alloy. Although relatively simple, it has enough features that could be used as a first test 
of the reliability of the technique. The main questions to be answered in this work relate to the 
existence of low temperature ordered structures for specific concentrations, for example, the ability 
to distinguish between rather similar phases for equiatomic alloys (CuAu I and CuAu II, the latter 
characterized by an antiphase boundary separating two identical phases). 


Engineering materials of necessity have various defects that result from annealing and pro- 
cessing procedures. In general, the ideal thermodynamic equilibrium is rarely achieved in the 
laboratory process of making the actual alloy, particularly in the case of large scale structures. It 
is of interest, therefore, if these procedures can be simulated. In this paper we have used Monte 
Carlo methods [4] to simulate this processing using the BFS method [1] for system energetics. We 
accomplish this end by starting from a random distribution at high temperature, then examining 
the structure by freezing the solid at various temperatures, and examining the micro-structure The 
organization of the paper is as follows: first we discuss the Monte Carlo procedure, then give a brief 
description of the BFS method and finally examine the defects appearing as a result of quenching 

thC "rte Monte Carlo procedure employed is a variant of that used to simulate the Ising system. 
The computational cell used in all simulations consists of 1008 atoms arranged on a face-centered 
cubic lattice. Boundary effects are minimized through the use of periodic boundary conditions. Each 
lattice site is characterized by the species of the atom (i.e. Cu or Au) occupying the site. The imUaJ 
state of the computational cell is a random alloy having a specified composition. The cell s lattice 
constant is assumed to depend only on composition and its variation with temperature is ignored. 
The lattice constant is determined from static zero-temperature BFS calculations (to be discussed 

^ A sequence of temperatures (the "cascade”) is chosen, in all cases, where, the system is allowed 
to equilibrate sequentially at each of the temperatures (simulating the "slow cooling of the actual 
alloy starting from a disordered solid solution) via a procedure where pairs of atoms of opposite 
species are chosen randomly, and the atomic species of each of the two atoms is reversed. The 
reversal is accepted or rejected using the standard Metropolis criterion. Because the lattice constant 
and all other properties are assumed to be independent of temperature, the temperature enters into 
the simulation only through the Metropolis criterion [4]. After the system has achieved equilibrium 
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(based on the total energy), the simulation is allowed to proceed further while various properties of 
the system are computed and averaged. These include the average energy of the cell, the specific 
heat, and bond correlations. 

While these simulations do not attempt to mimic the detailed dynamics of the equilibration 
process, they do offer a qualitative view of the effects of rapid versus slow cooling of the system. 
The temperature treatment (that is, the sizes of the steps between the various temperatures in the 
cascade) is of critical importance in determining the final state of the system. Slow cooling results in 
a highly-ordered low-temperature state, while rapid cooling results in an alloy with grain structure, 
in which each grain has essentially the ordering of the thermodynamic ground state. 

The current simulation procedure does not allow for relaxation of neighbors of atoms whose 
species have been reversed. In cases where the lattice constants of the two species are substan- 
tially different (which is the case in the Cu-Au system), neglecting such relaxation may result in 
a significant error in the system energetics, with the result that quantities such as phase trans- 
ition temperatures may not be accurately predicted. In order to eliminate this limitation and to 
more accurately predict transition temperatures, we have begun follow-up work on the Cu-Au sys- 
tem in which the lattice is relaxed via a molecular dynamics process, while the system is ordered 
concurrently using the Monte Carlo procedure described above. Because this new work is compu- 
tationally more intensive than that described in this paper, we are developing our computer code 
for a 16 -processor SP2 parallel processing computer using the PVM environment. 

The phase diagram of Cu-Au is shown in Fig. 1 [ 5 ]. A solid solution exists for the whole range 
of concentration, with a rather well defined order-disorder transition, that leads to simple LI2 (for 
CuAu3 and CU3AU) and Llo (for CuAu) fee ordered structures. Of special interest is the ordering 
found for equiatomic alloys. For a narrow range of temperatures, there is a long-period superlattice 
(CuAu II) where an antiphase boundary (APB) arises. While antiphase boundaries introduced by 
plastic deformation belong to the microstructure of the alloy, the former belongs to the phase. On 
either side of the APB, the alloy has the Llo structure, which for even lower temperatures becomes 
the only phase found for the Cu-Au alloy (CuAu I). The purpose of this work is to determine the 
ordered structures of Cu-Au as a function of composition by means of Monte Carlo simulations and 
to examine the influence of the temperature (heat) treatment in the simulation in relation to the 
type of ordering found. The energetics of the system are calculated with the BFS method for alloys, 
which has been successfully applied to the study of other properties of the Cu-Au system. In what 
follows we provide a brief description of the BFS method and some comments on the Monte Carlo 
calculations performed, followed by a discussion on the main features of the Cu-Au phase diagram 
obtained in this work. 

The BFS method is based on the idea that the energy of formation of an alloy is the superposition 
of individual contributions £; of non-equivalent atoms in the alloy [1]: 

£i=ef + s,-(£? -ef # ). (!) 

£i has two components: a strain energy e 5 , computed with equivalent crystal theory (ECT) [6], that 
accounts for the actual geometrical distribution of the atoms surrounding atom i, computed as if 
all its neighbors were of the same atomic species, and a chemical energy sP — £^° , which takes into 
account the fact that some of the neighbors of atom t may be of a different chemical species. For 
ef we interpret the chemical composition as a defect of an otherwise pure crystal. We represent 
this defect by ‘perturbing’ the electronic density in the overlap region between dissimilar atoms and 
locating them at equilibrium lattice sites of atom ». The ideas of ECT [6] are used to develop a 
procedure for the evaluation of the energy associated with this ‘defect’. To free the chemical energy 
of structural defect energy which should only be included in the strain energy, we reference ef to a 
similar contribution where no such perturbation is included (ef°). The coupling function <7,, which 
ensures the correct asymptotic behavior of the chemical energy, is defined as y, = e ~ a < , where af is 
a solution of ef = -E* c [l - (1 + af)exp(-af)] (see ref. [ 7 ]), and where ££ is the cohesive energy 
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for atom *. In the context of BFS, the terms ‘strain’ and ‘chemical’ represent quite different effects 
than the usually assigned meanings. For a dear understanding of this paper, we direct the reader to 
Ref. 1 where a detailed description of the calculation of the strain and chemical energy contributions 
is provided. Except for two parameters determined by fitting to experimental or theoretical alloy 
properties, the method relies on pure element properties. The BFS parameters used in this work, 
as well as the ECT parameters are listed in refs. 1 and 6 , respectively. The parameters used in this 
work have been previously used, with a great degree of success, in previous BFS applications to the 
Cu-Au system, induding surface segregation as well as multilayer surface relaxation studies [ 1 ]. 

There are some general observations regarding the results of the current study. Most simulations 
were performed following a ’cascade’ process: starting with a random configuration at very high 
temperature, the cell was stabilized for successive temperature steps, which could be loosely related 
to a slow cooling of the actual alloy starting from a disordered solid solution. The temperature 
treatment of the sample is essential in determining the final state. Slow cascade processes always 
resulted in highly ordered compounds at low temperatures, with very few antisite defects, whereas 
sudden cooling of the sample results in an alloy with grain structure, where each grain has basically 
the ordering pattern to be expected for the thermodynamic ground state. This process can be 
understood as if the quenching of the sample ’freezes’ ordered domains within the cell, whose seed 

is already present in the initial disordered state. 

Fig. 2 indicates the final states of ’cascade’ calculations for three specific concentrations: 25, 
50 and 75 at. % Au. As expected, at low temperatures these structures follow the Ll 2 (CuaAu and 
CuAu 3 ), Fig. 2 (a) and the Ll 0 (CuAu I) ordering Fig. 2(b). The coordination parameter <r, which 
indicates the probability that a Cu atom has a Au atom as a nearest-neighbor, should attain the 
values 0.3333 (CU 3 AU), 0.6667 (CuAu I) and 1.0000 (CuAu 3 ) if perfect order was achieved. In our 
calculations these values are 0.3327, 0.6630 and 0.9940, respectively. The final temperature for each 
one of these three alloys is T= 100 K.Slow cooling of certain disordered alloys to temperatures in the 
400 K range result in the CuAu II structure, with a clearly defined APB, as seen in Fig. 2 (c). 

Finally, as mentioned above, sudden cooling of a disordered alloy results in partially ordered 
alloys, as shown in Fig. 3. Close examination of the final cells shows the possibility of identifying 
domains within the computational cell which essentially follow the Llo ABAB stacking for the 
equiatomic case (Fig.3(b)) and the Llj ordering for alloys with 25 and 75 at. % Cu. To guide the 
eye. Fig. 3 shows a deformed computational cell for each concentration, where different (Fig.3( *).(«)) 
regions showing Llj ordering can be seen. The Ll 0 ordering can be seen in more than one direction. 
In order to highlight this fact, Fig. 4 shows a front view of the quenched CuAu alloy, where some 
domains are indicated as well as their orientation relative to each other. 

In spite of the microstructure of the alloys shown in Fig. 3 there is a high degree of ordering 
as measured by the probability of an atom A having a nearest-neighbor of species B. If A denotes 
a Cu atom, this quantities should have the values 1.0, 0.6667 and 0.3333 for CuAu 3 , CuAu, and 
Cu 3 Au, respectively. As mentioned above, a slow cooling process yields the values 0.9940, 0.6630 
and 0.3327. The corresponding probabilities for the samples shown in Fig. 3 are 0.9735, 0.6524 
and 0.3208, respectively, thus indicating that the domains of Llo *od Llj ordering locate themselves 
so as to optimize the expected coordination. The degree of order and the size and abundance of 
these domains is clearly a manifestation of the short range ordered that could exist in the high 
temperature alloy. 

Restrictions were imposed in these calculations although none of them, if lifted, are expected to 
lead to qualitatively different outcome. In all cases, the initial state for the Monte Carlo simulation 
was a random state at a very high temperature (5000 K), which is obviously well beyond the 
melting temperature. However, the atoms were constrained to locate themselves in equilibrium 
lattice sites. Moreover, no local atomic relaxations were allowed, therefore eliminating the possibility 
of detecting a tetragonal distortion. The computational cell was restricted to retain a constant 
volume throughout the calculation. The main effect of such a restriction could be the degree of 
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ordering achieved at each given temperature, but given the nature of the results presented here, it 
might be safe to assume that this limitation will not alter the essential nature of the final state, i.e. 
the type of ordering found. 

In conclusion, we have proposed the use of Monte Carlo techniques with the BFS method for 
energetics in order to simulate certain aspects of alloy processing for Cu-Au, namely, slow cooling 
and rapid quenching. This approach appears to reproduce some of the features of processing. Slow 
cooling to 100 K produces a single crystal for each alloy examined. Slow cooling to 400 K produces 
the CuAu II structure with antiphase boundaries as observed in the phase diagrams. Rapid cooling 
gives a granular micro-structure with grains having the equilibrium ground state structure of each 
alloy studied. This Monte Carlo procedure could then be a possible non-dynamic method for 
studying the effects of alloy processing. 
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Fig. 1: Phase diagram for the Cu-Au system [5]. 



Fig. 2 Final geometries for slow cooling 
to 100 K giving the correct ground state 
structures of (a) CuAu3 (LI2) (or CU3AU) 
(b) CuAu I (Ll 0 ) and (c) CuAu II. Solid 
spheres denote Cu atoms. 
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Fig. 3 Expanded cells for quench cooling to 400K: (a) CuAu 3 , (b) CuAu, and (c) Cu 3 Au showing 
resulting grain structure. Solid spheres denote Cu atoms. 
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Fig. 4 End view of Fig. 3(b) showing domains with ABAB structures in each direction following 
quench cooling to 400K. Solid spheres denote Cu atoms. 
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